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Abstract. A cyclic random walk is a random walk whose transition probabilities/rates 

can be written as a superposition of the empirical measures of a family of finite cycles. 

This identifies a convex set of models. We discuss the problem of characterization of cyclic 

»vj ■ random walks in some special cases showing that it is related to several remarkable and 

classical results. In particular we introduce the notion of balanced measure and show 

^^ ■ that a translation invariant random walk on Z'' is cyclic if and only if its transition 

^SJ ' probability is balanced. The characterization of the extremal elements is obtained using 

> I , the Caratheodory's Theorem of convex analysis. We then show that a random walk on a 

O f finite set is cyclic if and only if at every vertex the outgoing flux of the transition graph 

is equal to the ingoing flux. The extremal elements are characterized by the Birkhofi^- 

Von-Neumann Theorem. Finally we consider the discrete torus and discuss when the 

cyclic decomposition can be done using only homotopically trivial cycles or elementary 

cycles associated to edges and two dimensional faces. While in one dimension this is 

equivalent to require some geometric properties of a discrete vector field associated to 

^^ ■ the transition rates, in two dimension this is not the case. In particular we give a simple 

^ , characterization of the polyhedron of the rates admitting a cyclic decomposition with 

• ■ elementary cycles. The proof is based on a discrete Hodge decomposition, elementary 

1 -^ I homological algebra and the Helly's Theorem of convex analysis. Finally we discuss a 

j^ ■ natural discretization procedure of smooth divergence free continuous vector fields and 

an application to random walks in random environments. 
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Pr . 1. Introduction 

A discrete/continuous time cyclic random walk is a random walk whose transition prob- 
abilities/rates can be written as a superposition of empirical measures of a family of finite 
___^_ cycles. This definition identifies a convex set of models. The aim of this paper is to 

S^ ■ characterize this set in some special cases. This is important since once you known that a 

H , random walk is cyclic then several properties and results can be deduced. For example the 

following results have been obtained using some kind of cyclic decomposition: in [9] and 



[T2j a central limit theorem and in ^ a quenched central limit theorem for random walks 
in random environments, in [T7j some regularity properties of the diffusion coefficient for 
a mean zero exclusion process, in [18] a central limit theorem for a tagged particle of 
a mean zero exclusion process, in [13] some bounds on the transition probabilities. In 
different papers some slightly different cyclic decompositions are used. In any case all of 
them correspond to define, starting from the Markov model, a weighted oriented graph for 
which a cyclic decomposition holds. Our results can be interpreted as results concerning 
cyclic decompositions of weighted oriented graphs. In this sense they can be applied in 
any case independently from the special rule associating a oriented weighted graph to a 
Markov model. There is an huge amount of results concerning cyclic decompositions of 
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graphs, see for example [T], [3] and [6]. The majority of them are of combinatorial type. 
Our results have instead a geometric flavour and are the following. 

We consider translational invariant Markov chains on the lattice Z . A model is iden- 
tified by a probability measure on Z which determines the weight to be associated to 
the different jumps. We introduce the notion of balanced measure. An Hyperplane on 
M containing the origin determines two half spaces. We can then compute the average 
distance from the hyperplane for the measure restricted to the two half spaces. If the two 
averaged distances are equal (possibly also +00) and if this happens for any hyperplane 
then we say that the measure is balanced. Clearly any mean zero measure is balanced. 
A translational invariant Markov chain is cyclic if and only if the corresponding measure 
determining the distribution of the jumps is balanced. A characterization of the extremal 
elements is obtained using the classic Caratheodory's Theorem of convex analysis whose 
statement is recalled in section [31 The result for mean zero measures could be deduced by 
the results in pT] and [21J. Here we give an independent proof and extend the validity of 
a cyclic decomposition to the class of balanced measures. This result is not a special case 
of the Choquet-Bishop-de Leeuw theorem [TS] since compactness is missing. See [TU] and 
[2U] for a general discussion of problems of this type. 

On a finite set a Markov model is determined by an oriented weighted graph. An 
oriented weighted graph is said to be balanced if for any vertex the ingoing weight is equal 
to the outgoing weight. We show that a Markov model on a finite set is cyclic if and 
only if the corresponding oriented weighted graph is balanced. In the case of discrete time 
Markov chains, a characterization of the extremal elements is obtained using the classic 
Birkhoff-Von Neumann Theorem whose statement is recalled in section [3l The proofs of 
these results are elementary but are useful in the discussion of the next issue. In [2j we will 
discuss and apply a generalization of this result valid for infinite graphs with a condition 
of vanishing flux towards infinity. 

On a finite graph we can consider a restricted class of cycles taking into account some 
topological obstructions. More precisely we can consider homotopically trivial cycles or 
elementary cycles associated to one and two dimensional faces of some cellular decomposi- 
tion. We discuss this problem in the cases of the one dimensional and the two dimensional 
discrete torii. While in one dimension the validity of restricted cyclic decomposition of 
this type is equivalent to require some geometric properties of a discrete vector field con- 
structed starting from the weights, this is not the case in two dimensions. In particular 
we give a simple characterization of the models that admit a cyclic decomposition using 
elementary cycles. Similar computations can be done also in higher dimensions. In these 
cases it is necessary a more detailed discussion of the homological structure and a more 
involved geometric construction. We will discuss this generalization in a separated paper 
[7] together with a discussion of the case of infinite grids.. 

We discuss also some applications. First we introduce a simple and natural discretiza- 
tion procedure for a continuous smooth divergence free vector field. The result is a discrete 
divergence free vector field on the lattice. Then, using the previous results, we give a con- 
dition to obtain a weighted oriented graph admitting an elementary cyclic decomposition. 
Finally we apply the above construction to deduce a quenched Central Limit Theorem for 
a random walk in random environments using the results in [5]. 

The paper is organized as follows. In section [2] we fix notation and state our main 
results. The section is subdivided into four subsections according to the three different 
frameworks that we discuss plus a subsection for applications. In section [3] we recall some 
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classic results and prove some original results that are useful in the proofs of our main 
results. In section H] we collect the proofs of our main results in the case of a translational 
invariant Markov chain on the lattice Z . In section [5] we collect the proofs of our main 
results in the case of a Markov model on a finite set. In section [6] we collect the proofs 
of our main results concerning cyclic decompositions using restricted classes of cycles. In 
section [7] we discuss the applications. 

2. Notation and main results 

In this section we state all the main results of the paper and fix notation. 

2.1. Cyclic random walks on JJ^. An element x G Z'^ has coordinates (xi, . . . ,Xrf). We 

call e*-*^ the versors of the canonical basis of M . This means that we have e^* := (Jjj, 

where b is the Kronecker delta. We denote by x • y := Yli=\ ^iVi ^^^ Euclidean scalar 
product among x,y G W^. With |x| := ^J^^ we denote the Euclidean norm of x. 

We denote by A^ ^ ^ the set of positive measures on U^ with total mass less or equal to 
1 endowed with the topology of weak convergence. Since in all the cases we will consider 
a tightness condition will be always satisfied, the weak convergence will be equivalent to 
the pointwise convergence. In A^ ^ ^ there is a natural partial order structure. Given 
p, (7 € 7W ^ ^ then p ^ q when p{x) ^ q{x) for any x G IJ^. The support of p G 7W ^ ^ is 
defined as 

S{p) := ixel'^ : p{x) > 0| . 

By A^^ (1 M.^^ we denote the subset of probability measures and by A^o ^ A^ ^ ^ the 
subset of mean zero measures. This is the set of measures p such that 

p{x)\xi\ < +00, Vz, (2.1) 

and moreover Ylx^l'^Pi^)^ ~ ^• 

Given A C M*^ we denote by co{A) its convex hull and by aff{A) its affine hull. These 
are defined as 



E 



n n 



co{A) := <x eM.'^ : X = '^ aix"; x* G A, Oj ^ 0, ^ a^ = 1, n G N 



i=x i=x 

n n 



aff{A) := <^ X G M'^ : x = ^ QiX*; x* G A, a^ G M, ^ a, = 1, n G N ^ . 
Let z^, z^, . . . , z^~^ be distinct elements of Z and consider also z" := z^. Define 

as the finite cycle on Z*^ that, starting from z^, visits sequentially all the elements z* in 
increasing order with respect to the index i. We say that such a cycle C has cardinality 
\C\ = n. Given a finite cycle C we define its empirical measure p G A1^ as 

1 ""' 

where 6x is the delta measure concentrated at x. The cycles of cardinality 1 are of the 
type C = {z^,z^) and the corresponding empirical measure is Sq. Empirical measures of 
finite cycles are called purely cyclic measures. It is clear that p defined in ()2.2p depends 
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only on the values of the displacement vectors w^ := z^ — z^ ^ and not on their relative 
order. Note also that for any cycle C we have 

ic^l \c\ 



=1 j=i 



As a consequence we have 



1 ""' 

so that every purely cyclic measure has mean zero. 

Consider a time homogeneous, discrete time random walk {X„}^gpj on Z with transla- 
tion invariant transition probabilities determined by p G TW ^ by the following identification 

P(X„+i = y\Xn = x) := p{y - x) . (2.3) 

A translation invariant discrete time random walk on Z, is called purely cyclic if the 
measure p in (|2.3p is purely cyclic. This means that there exists a finite cycle C such that 

P = P . 

As far as p is concerned, we can naturally introduce an equivalence relation ~ between 
cycles. Let C be a cycle and let {w^, . . . , w''-''} be the corresponding set of displacement 
vectors. Note that in the set of displacement vectors of a cycle, a vector w may appear 
more than once. For this reason we will also denote the set of displacement vectors as 

{{w^,ni),... ,{w'',nk)}, 

to indicate that the vector w^ appears nj G N times. Consider now another cycle C and 
let {v^, . . . ,v''-^ 1} be the corresponding set of displacement vectors. We say that C ^ C 
if 

{w\...,w\''\] = {v\...,v\^'\]. 

The equivalence class to which the cycle C belongs is denoted by \C] and is determined 
by the set of the displacement vectors of any representant C. We can then naturally use 
the following identification 

[C]^{{w\n^),...,{w\nk)}. (2.4) 

We denote by C the countable set of equivalence classes of cycles, endowed with the 
discrete topology. Note that if {{w^,ni), . . . , {w^ ,nk)} G C then 

fc 
^ni'u;*=0. (2.5) 

i=l 

Conversely if (j2.5p holds then there exist natural numbers n[ such that n^ ^ n^ and more- 
over {{w^ ,n'i), . . . , {w^,n'f^)} G C. We cannot use directly the numbers Ui since according 
to our definition a cycle is self-avoiding. Observe that if C ~ C" then p = p and 
we can therefore use the notation p^^^ to denote p*^ where C is any representant of [C]. 
The converse is not necessarily true. Indeed let / G N, / ^ 2, and consider any cycle C 
satisfying ()2.4p and any cycle C" with 

[C'\ = {{w\lni),...,{w\lnk)]. 

Then clearly it holds [C] 7^ [C] but nevertheless p = p . 
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Definition 2.1. A probability measure p on TL is called cyclic if there exists a probability 
measure p on C such that p = pP where p^ is defined as 

pP := Y, p([C])p[^l . (2.6) 

[C]ec 

A translation invariant discrete time random walk on TL'^ is called cyclic if the measure p 
in ()2.3p is cyclic. 



The meaning of (12 .Sp is that for any x € Z we have 



pp{x)= ^p([C])p[^](x), (2.7) 

[C]ec 



i.e. we require just pointwise convergence. By monotone convergence Theorem ([27 
imphes 

pPiA)= Y p{[C])p^^\A) , \/ACZ^. 

[C]GC 

Equivalently we have a probabihty measure on C x Z"^ that gives weight p{[C])p^'^' {z) to 
the pair ([C],2;) and p^ is its Z'^ marginal. Definition (|2.6p can also be interpreted as 
follows. For any sequence C„ C C such that \Cn\ < +oo, C„ C C„+i and C = UnC„ it holds 

p''= hin Yl /'([C])P"'^ (2.8) 

[c]eCn 

where the limit is in the weak sense. This follows immediately by pointwise convergence 
since a tightness condition is clearly verified. 

If the measure pP is integrable, i.e. if (12. ip holds with pP instead of p, then by Fubini 
Theorem we get 

Yp'i^)^= E ( E p(P]y^'(^) ) ^= E ^([^]) f E ^p^'^Hx)] =0. 

xezd xezd \[C]£C J [C]ec \x€Zd J 

This means that every measure pf that is integrable has mean zero. 

Given H an hyperplane in R*^ we call H^ and H^ the two closed half-spaces determined 
by H. Given A,B CR<i we define 

d{A,B):= inf |x — y| . 

{xeA,yeB} 

Definition 2.2. A measure p ^ Ai^^ is called balanced if for any hyperplane H B we 
have 

E p(.x)d{x,H)= E p(.x)d{x,H). (2.9) 

X&H+ x£H- 

We call B <^ Ai^ the set of balanced probability measures. 

The main result of this subsection is the following. 
Theorem 2.3. An element of Ai^ is cyclic if and only if it is balanced. 



We underline that in ()2.9p both sides can also be infinite. Note also that MoHAi Q B. 
Indeed let n be a versor normal to H then we have H~^ = {x : x ■ n ^ 0} and H^ = 
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Figure 1. An example in 7? of a measure p that satisfies condition ()2.9p 
with respect to the 2 hyper-planes orthogonal to e^^' and e^"^' but it is not 
balanced and consequently also not cyclic. With • we denote elements of 
5(p), the dashed region represents co{S{p)). Condition ()2.9p is not satisfied 
for the hyperplane H. 



{x : X • n ^ 0}. Given x £ H~^ it holds d{x, H) = x-n and likewise given x € H it holds 
d{x, H) = —X ■ n. If p G Mq n A^"*^ then we have 

0= y^ p{x) {x ■ n) = y^ p{x)d{x,H) — y, p{x)d{x,H) . 

Indeed the same computation shows that when p is integrable then it is balanced if and 
only if it belongs to A^o H A^^. This corresponds to verify condition ()2.9p just for the d 
hyper-planes orthogonal to the vectors e^'^' of the canonical basis. The validity of ()2.9p 
for the others hyperplanes will follow. When p is not integrable the balancing condition 
(|2.9p has to be checked for any hyperplane H, since the validity of ()2.9p with respect 
to the hyper-planes orthogonal to the vectors of the canonical basis does not imply the 
validity of (|2.9p for any H. Consider for example the measure p defined by p{2i, —i) = 
p{—i, 2i) = -^ , i ^ 1 and giving weight zero to all the remaining elements of Z^ (see 
figured]). This measure satisfies ()2.9p with respect to the 2 hyper-planes orthogonal to the 
vectors e^^^e^"^' but it does not satisfy the balancing condition for example with respect 
to the hyperplane H in figure [TJ This measure is not cyclic. Indeed it is not possible to 
generate a cycle using vectors in S{p) since co{S{p)). 

A simple example of a balanced not integrable probability measure on I? is p(x) = — ^ 

Another example is given by a probability 



where c is a suitable normalization constant. 
^^ such that 



measure on 



+00 . 



^ p{{xi,0))\xi\ = Y^ p((a;i,0))|a;i| = ^ p{{0,X2))\x2\ = ^ p((0,X2))|x2 
xi > xi ^ X2 > X2 ^ 

Then the above conditions guarantee that the measure is balanced and consequently cyclic 
whatever are the values assumed by the measure outside from the coordinate axis. 

We discuss now in detail the one dimensional case. In one dimension a balanced prob- 
ability measure p either belongs to A^o H 7W^ or satisfies 



y xp{x) = y^ xp{—x) = +00 . 



(2.10) 



x=l 



x=l 
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We now sketch an iterative procedure whose generahzation is at the core of the proof of 
Theorem 1 2 . 3 1 and will be discussed more in detail during its proof. This procedure generates 
a cyclic decomposition of p satisfying ()2.10p . We remark that this case corresponds to a 
not integrable balanced measure. Given a measure p € A4^^ satisfying ()2.10p we define 



X' 



minjx ^ 1 : p (x) > O} , 



x_ := max {x ^ — 1 : p (x) > O} . 

We need to distinguish two cases 

p'(x'^)x'^ + p\x^_)x^_ ^ , (a) , 

p\x^^)x^j^ + p\x^_)x^_ < , (6) . 

We consider the equivalence class of cycles 

where n^ and n_ are determined requiring that —f- is an irreducible fraction and moreover 
-r = — p. The corresponding weight will be 

— 'iT: — (")' 

v\x\){n\W_) ,^ 

If we define p'+^ := p' - p([C'])p['^'l then p'+^ G A4 ^ ^ and stih satisfies (|2.10p so that 
we can iterate the procedure starting from p^ := p — p(0)(5o. By construction we have 
I {^S{p^'^'^) n [x'_,x',_]} I ^ 1 and condition (j2.10p guarantees that lim/^_|_oo 2;- = —00 and 
lim/^+00 x\ = +00. As a consequence we get 

+00 

1=1 
that is a cyclic decomposition like (|2.6p . 

We call irreducible an equivalence class [C] = {w^ , . . . w^^^} such that for any non empty 
{v^ , . . . v^} C {w^ ^ . . . w''-^'}, where the inclusion is strict, it holds Y2i=i ^* / 0- 

In lemma 14.11 stated and proved in section HI we will characterize C* a subset of the 
irreducible equivalence classes of cycles. The second main result of this section is the 
following. 

Theorem 2.4. The set of extremal elements of B is {p }ic]pc*' Moreover for any 
element p & B there exists a probability measure p on C* such that 

P= Yl P(iC])P^''^- (2.11) 

[C]6C* 

We give also a different characterization of cyclic measures. We denote by A^o the 
monotone non-decreasing weak closure of A^o- This means that p € Mo if there exists a 
sequence pn G M.o such that pn ^ Pn-i and p = lim„^+oop„ in the weak sense. Then we 
have the following alternative characterization. 
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Theorem 2.5. 

The monotone requirement is necessary since it is easy to see that Mq n A^"*^ = Ai^, 
where A^o is the weak closure of TWq- 

In this subsection we restricted our discussion to probabihty measures but the arguments 
could be adapted to positive measures having finite mass. 

2.2. Cyclic random walks on a finite set. Both a discrete time and a continuous 
time Markov chain on a finite or countable set V can be described in terms of a oriented 
weighted graph (V, E, r). Here V is the set of vertices, E C V x V is the set of oriented 
edges and r : E ^ 1R+ is the weight function. For simplicity here and hereafter we consider 
the case (x, x) ^ E for any x £ V, the general case can be easily handled. The set of edges 
E is determined by the weight function by the requirement that {x,y) G -E if and only if 
r(x, y) > 0. For this reason we can also denote a weighted graph simply by the pair {V, r). 
The set of edges having positive weight is denoted by E{r). There is a natural partial order 
among weighted graphs having the same set of vertices. We say that (V, r) :< {V, r') if 
r(x, y) ^ r'{x, y) for any (x, y) G E{r). The set of all non negative weights r : V xV ^ M+ 
is denoted by W. 

In the case of a discrete time random walk we fix the weights in such a way that 

r(x, y) := F{Xn+i = y\Xn = x) . 

In this case clearly the normalization condition 

Y^ rix,y) = l, yx£V, (2.12) 

{y£V:{x,y)&E{r)} 

has to be satisfied. In the case of continuous time Markov chains we identify the weight 
r{x,y) with the rate of jump from x to y. 

A cycle C, on an oriented graph (y,E), is a finite sequence C := {x^,x^, . . . ,x"~^,x'^) 
of distinct elements x* G y such that for any z = 0, . . . , n — 1 we have (x*, x*"*"^) G E where 
the sum in the indices is modulo n. We say that an edge (x, y) belongs to the cycle C and 
write (x,y) G C if there exists an i such that (x,y) = (x*,x*^^). In this framework it is 
natural to identify two cycles if they contain the same set of edges. More precisely given 
two cycles C = (x°,x^, . . . ,x"'~-^,x°) and C = {y^,y^, . . . ,y"'~^,y°), we say that C ^ C 
if there exists a j G {0, 1, . . . , n — 1} such that y* = x*"'"-' ,i = 0, l,...,n — 1, where the sum 
in the indices is modulo n. As before we call [C] the equivalence class of cycles having C 
as a representant. We call also C the set of equivalence classes of cycles. 

We associate to the cycle C the weighted graph, that we call a purely cyclic weighted 
graph, determined by the weight function 

.c/^ „N ._ / 1 if (a;,y) gC, 



^ (-'^)-io othe;;;se. ' (2-13) 

Note that if C '^ C then r = r so that we can use the notation A ' without ambiguities. 
Moreover a purely cyclic weighted graph uniquely determine an element of C 

Definition 2.6. We call a weighted graph {V,r) cyclic if there exists a positive measure 
p on C such that 

r = Yl /'([C])rf^l . (2.14) 

[C]eC 
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Definition 2.7. We will say that the weight r is balanced at x £ V if it satisfies the 
condition 



E 



r{y,x)= ^ r{x,y). (2.15) 

{y:{y,x)eE(r)} {y:(x,y)eE(r)} 

The weight r is balanced if ()2.15p holds at every x £ V. 

We use the term balanced with two different meanings, one determined by definition 
12.21 and the other one determined by definition 12.71 There is no risk of confusion since one 
refers to elements of A^ ^ ^ and the other one to elements of W. The main result of this 
subsection is the following. 

Theorem 2.8. A finite weighted graph {V, r) is cyclic if and only if it is balanced. 

Remark 1. Note that all the results of this paper can be formulated in terms of cyclic 
decompositions of weighted graphs without any reference to Markov models. A natural way 
to associate a weighted graph to a Markov model is the one discussed at the beginning of 
this subsection but it is clearly not the only one. For example given a positive measure vr 
on V and a random walk with jump rates c(x, y) we can associate to this pair a weighted 
graph with weights 

r{x,y) ■.= 7r{x)c{x,y). (2.16) 

The definition of centered random walks in [13] is based for example on the properties of 
a weighted graph defined in this way. All the results of the present paper can be clearly 
applied also in this case. Note that the validity of the balancing condition of Definition \2. 1\ 
for the rates in (|2.16|) is equivalent to require that vr is invariant for the rates c. 

2.3. Cyclic random walks on a finite graph with topology. We consider the d 
dimensional continuous torus T = M /Z of side length 1. Any element of the torus is an 
equivalence class that we identify with any of its representant. On the continuous torus we 
embed a d— dimensional discrete torus with mesh -^. This is a graph whose set of vertices 
is 

VN:=lT'^3x = {xi,...,Xd) : xi, . . . , x^ = 0, ^, . . . , ^^^| . (2.17) 

The set of oriented edges E^- contains all the pairs {x,y), with x,y £ Vn and such that 
d{x, y) ^ 1/-^, where the distance is on T . For simplicity we consider a cubic grid but 
all the arguments can be repeated also in the case of a grid with A'^i x A'^2 x • • • x A^^ sites. 

We define E'j^ C Ej^ the collection of oriented edges {x,y) G Ejy such that there exist 
X, y G M in the same equivalence classes respectively of x, y and such that y = x + e^'^' /N 
for some i. Note in particular that for example ((— ^, X2, • • • , x^), (0, X2, . . . , x^)) G E'^. 
This is the exact meaning that we give to expressions like (x, y) G -Eat with y = x^e^"^' /N. 

The embedding of the discrete torus on the continuous one defines in a natural way a 
cellular decomposition of T'^. Every cell of maximal dimension is a d dimensional hyper- 
cube of side length -j^ having vertices in V]\j. Every 2 dimensional cell of the cellular 
decomposition is a square having vertices 

|x, X + e^'^/N, X + e^J'ViV, x + e'^^/N + e^^'V^} , 

for some x G Vat and some i ^ j. Every 2 dimensional cell can be oriented in 2 possible 
ways. An orientation is a choice among the two possible elementary cycles 

(x, X + e^^^/N, X + e«/^ + e^^^/N, x + e^^'^/N, x) , , . 

{x,x + e^i^ /N ,x + e^^N + e^^'^ /N ,x + e^''^ /N ,x) , ^ ' ^ 
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(a) (6) (c) 

Figure 2. Elementary cycles associated to one dimensional cells (a) and 
two dimensional cells (6) and (c), respectively clockwise and anticlockwise 
oriented. 

associated to the 2 dimensional cell (see figure [2] cases (b) and (c)). We call -F/v the set 
of oriented two dimensional cells of the cellular decomposition. Conventionally the two 
possible orientations are called clockwise and anticlockwise. If / G -Fat we call f^ the 
element of F/v associated to the same 2 dimensional cell but corresponding to the opposite 
orientation. We call F'j^ C -F/v the subset of anticlockwise oriented 2 dimensional cells, 
this means the cells with associated the orientation 

(x, X + e(™>^^*'J'»/iV, X + e^'^/N + e^^^/N, x + e^'^^^^*'^'"/^^, x 

The discrete torus is a finite graph so that cycles and equivalence classes of cycles are 
defined as in subsection 12.21 In particular we still call C the set of equivalence classes of 
cycles. Given a cycle C = {z^, z^, . . . , z"'~^, z") on the discrete torus we associate to it the 
displacement vectors w^, . . .w^ defined by 

% '.:'.-i!^' (2-19) 

TV n z — z ^ . 

Given a cycle C = {z^, z^, . . . , z''^~^, z^) on the discrete torus we can naturally associate 
to it also a continuous closed curve on T*^. This is the projection on T'^ of a continuous 
curve on M'^ whose parametrization {x(s)}^grg ^i is defined by the positions x(0) = z^ and 



s e 



1 i 



n n 



Since the projection is a closed curve on T then x(l) — x(0) is a vector with integer 
coefficients. A continuous closed curve on T is called homotopically trivial if it can 
be continuously deformed to a single point. More precisely a continuous closed curve 
{x(s)}ggrQ ;^i is homotopically trivial if there exists a continuous map X : [0, 1] x [0, 1] — >• T'' 

such that X{s,0) = x{s), X{s, 1) = x* for any s G [0, 1], where x* G T'^ is a fixed element 
and X{l,t) = X{0,t) for any t. The discrete cycle C will be called homotopically trivial 
if the associated continuous curve on the torus is homotopically trivial. It is easy to see 
that this condition is equivalent to require that 

n 

^u;* = 0, (2.20) 

1=1 

where w^ are the displacement vectors of the cycle. By C* C C we denote the subset of 
equivalence classes associated to homotopically trivial cycles. 
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Figure 3. Elements of Fj^ of the two dimensional torus. Any pair is 
oriented in agreement. Opposite sides of the square are identified. 



We define also a subset C^ C C* of equivalence classes of elementary cycles. Elementary 
cycles are naturally associated to one and two dimensional cells of the cellular decompo- 
sition. Given x,y £ Vjy such that d{x, y) = 1/N we call [Cj^, ,^}] the equivalence class 
containing the elementary cycle C^^yj := {x,y,x) (see Figure [2] case (a)). 

As already discussed to every element / E Fj\[ corresponds a cycle of the type (j2.18p . We 
call [Cf] the corresponding equivalence class. The set of equivalence classes of elementary 
cycles C^ is defined as 



C' 



{^{^•,S/}]}( 



u {[C/]} 



/eFjs 



Given (x, y) G Ej\f and / G F/v we write {x, y) G / if the oriented edge {x, y) belongs to 
any representant of [Cf]. A pair of elements of Fjy associated to adjacent cells is said to 
be oriented in agreement if no elements of E^ belong to both. On the two dimensional 
torus every pair of elements of F^ is oriented in agreement (see figure [3|). This is possible 
since the two dimensional torus is an orientable surface. 

Given a continuous time random walk on the discrete torus its transition rates identify 
a weight r G W. The discrete torus is a finite graph and a decomposition like (j2.14p exists 
if and only if the balancing condition of Definition 12.71 holds. We investigate under which 
conditions decompositions like 



[C]GC* 



or 






p{[C])r 



[c] 



(2.21) 



(2.22) 



hold. Recall that p is a positive measure and A ' is the element of W associated to any 
cycle in [C] like in ()2.13p . Let us call i?* C VF and R^ (IW the sets of weights for which 
respectively a decomposition like (|2.2ip and like (|2.22p holds. Clearly it holds W C R* . 

Before proceeding we need to introduce some notions and terminology. We refer to 
[l]j [ID]) |14] for more details. We call A'^ the |y/v|— dimensional vector space of the real 
functions / : I^at — ^ M. 

We call a function (p : E]\f — t- M a discrete vector field if it satisfies the condition 



(t>{^,y) = -(t}{y,x) 



V(x,y) G En 



(2.23) 
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We call A^ the l-E^I— dimensional vector space of discrete vector fields. We endow A of 
the scalar product (•,•)! defined as 

(0, </.')! := Y. <Pix,y)<P'ix,y)- (2-24) 

Due to ()2.23p the value of ()2.24p does not depend on the specific choice of E'j^. On W we 
define a projection operator that associates to a weight r a discrete vector field cp"^ G A^ 
defined as 

(p''{x,y) ■■= r{x,y) -r{y,x). 
Conversely to any (/> £ A^ we canonically associate a weight r'l' £ W defined as 

r'^(x,y) := [(/>(x,y)]+, 

where [•]+ denotes the positive part of a real number 

a if a ^ , 
otherwise . 

Note that min |r'^(x,y),r'^(y,x)} = and r"^ G -R(0), where 

R{(j)) ■= {r € W : c/)^ = (/)} . 

We stress that R{0) coincides with the set of symmetric rates, i.e. such that r{x,y) = 
r{y, x). Given a positive weight r we can decompose it into the sum of two positive weights 
as 

r = s + r'l'\ (2.25) 

Clearly s G ^(0) and it is called the symmetric part of r. We can deduce 

R{(l)) = |r : r = s + r'^ , s£ R{0)\ = r't' + R{0) . (2.26) 

A simple consequence of (j2.26p is that there exists a minimal element, with respect to the 
partial order <, on R{(j)) that is r^. 

We call a two chain a map ip : -F/v — ^ ^ that satisfies the additional property 

Hf) = -Hn, V/GF;v. 

We call A^ the |i^J^|— dimensional vector space of the 2-chains. 

We define the co-boundary operator 5 : A* — )• A*+^, i = 0, 1 as follows. Given a function 
/ G A° we define 6f G A^ as 

6f{x, y) := f{y) - f{x) , (x, y) e En ■ 

A discrete vector field obtained in this way is called a gradient. Given </) G A"*^ we define 
6(j) G A^ as 

mf)-= E 't>{x,y), feFN. 

{ix,y)GEN:(x,y)ef} 

The value 64>{f) is called the value of the circulation of the vector field (p around /. 

We define also a boundary operator d : A* — )• A*~^, i = 1,2. Given ^ G A^ we define 
dijj G A^ as 

#(x,y):= Yl ^(/)' ix,y)&EN. (2.27) 

{/6Fjv : {x,y)€f} 
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Given (/> £ A^ we define d(j) S A'^ as 

d4>{x):= Y. H^^y)^ xeVN. (2.28) 

{vGVn ■■ {x,y)eEN} 

The r.h.s of (j2.28p is called the discrete divergence at x of (j). 

It is easy to verify according to our definitions that for any -0 G A^ it holds 

d(#) = 0. (2.29) 

The one dimensional case is elementary and will be considered separately. Here we 
restrict to dimensions d ^ 2. We denote by dA^ C A-*^ the set of vector fields (j) for which 
there exists a ip £ A'^ with 4> = dtp. A characterization of dA^ will be discussed in terms of 
a discrete Hodge decomposition in section [3l It can be summarized as follows. A discrete 
vector field (j) belongs to dA^ if and only if d(j){x) = for any x and 

Y^ (p(x,x + e'-''>/N] =0, i = l,...,d. 

x£Vn 

We have the following result. 

Lemma 2.9. A necessary condition to have r € R* is that (p^ £ dh? . The same holds for 

We have also the following Lemma where the metrics in W is the Euclidean one. 

Lemma 2.10. The sets W and R* are closed sets in W . 

As we will see, an immediate consequence of the above Lemma is the following Corollary. 
It says that (p"^ G dh? is not a sufficient condition to have r £ R* (and consequently also 
r G R^). We add a condition of positivity since in this way we get a stronger statement. 

Corollary 2.11. There exists a weight r £ W such that (ff E dh? , r{x,y) > for any 
{x,y) G En hut nevertheless r ^ R* . 

2.3.1. One dimensional torus. In the one dimensional case there are no two dimensional 
faces and 

so that \C*\ = N. Moreover we have \C\ = N -\- 2. To the previous equivalence classes we 
need to add the two equivalence classes having as representants the cycles 

C+:={0,1/N,2/N,...,{N-1)/N,0) 

and 

C. := (0, (N - l)/iV, (N - 2)/N, ..., l/N, 0) . 

The following Theorem says that the validity of decompositions of the type ()2.2ip and 
()2.22p is equivalent to a geometric property of the vector field (p^' . More precisely ()2.2ip 
and ()2.22p hold if and only if (j)^ is zero. 

Theorem 2.12. It holds a decomposition of the type ()2.2ip if and only if(p'' = 0. Moreover 
this happens if and only if the rates r are reversible with respect to the uniform measure. 

Let 

m := min r{x,y) . 

(x,y)eEN 

The following Theorem says that also the validity of a decomposition like ()2.14p is still 
equivalent to a geometric property of the vector field cp"^ . More precisely ()2.14p holds if 
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and only if (jf is divergence free. Moreover the model is so simple that we can characterize 
completely the associated measures p. 

Theorem 2.13. A decomposition like (j2.14p holds if and only if the discrete vector field 
(ff' is constant or equivalently has zero divergence dcff = 0. Moreover if cjf = c then all the 
measures on C for which ()2.14p holds are parameterized by the real parameter a G [0, m] 
as 

P(.[C{x,y}]) = min {r{x, y),r{y, x)} - a , 

p{[C+]) = [c]+ + a, (2.30) 

p{[C.]) = [-cU + a. 

By a constant discrete vector field (p = c we mean (j){x, y) = c for any (x, y) G E'j^. 

2.3.2. Two dimensional torus. Every (x, y) G E'^ belongs to only two elements of -F/v and 
moreover one is clockwise oriented and the other one is anticlockwise oriented. We call /+ 
the element of F'j^ such that (x,y) G /+ and /__ the element of F'j^ such that {y,x) G /_. 
Let (f) G dl^ and consider any ip ^ I^ such that dip = 4>. As will be explained more in 
detail in section [3l any other ip' such that dip' = <j) differs from ip by an additive constant, 
this means that there exists a real constant c such that 

V''(/) = V'(/) + c, V/GF^. 

The main result of this subsection, Theorem 12.141 below, will not depend on the specific 
choice of c. To every (x, y) G E'j^ we associate the closed interval /(x, y) of R defined as 

I{x,y) = [ii{x,y),i2{x,y)] := [min{V'(/_), V'(/+)} ,max{V'(/-),^(/+)}] . 

Given /i, . . . ,/„ a collection of closed intervals of M we call V{Ii, . . . ,In) ^ IK" the 
closed unbounded convex polyhedron defined as follows. An element s = (si, . . . , s„) G M" 
belongs to P(/i, . . . , /„) if and only if the inequalities 

Si + Sj ^ d{Ii,Ij) , i,j = l,...,n, (2.31) 

are satisfied. Note that in ()2.3ip we are considering also the cases i = j that imply Sj ^ 0. 
See Figure m for a simple example in dimension n = 2. Two important properties of this 
polyhedron are the following. The first one is that if s G 'P{Ii, . . . ,In) and s' ^ s then 
s' G V{Ii, . . . , In)- The second one is that G Vih, ■ ■ ■ , In) if and only if d{Ii,Ij) = for 
any i,j. By the Helly's Theorem recalled in section [3] this is equivalent to Dili ^ 0. In 
this case clearly P(/i, . ..,/„) = (M+)'^. 

The main result of this subsection is the following. It says that, differently from the one 
dimensional case, the validity of a decomposition like ()2.22p is not simply equivalent to 
require a geometric property of the vector field (j)"^ , but involves instead also the symmetric 
part s. 

Theorem 2.14. Let r such that (ff = (j) ^ dh? . We have r ^ R^ if and only if 

sGP({/(x,y)}(,^j^)g^j^) . (2.32) 

In (|2.32p s is the symmetric part of r (see ()2.25p ) and the intervals are constructed using 
any ^ G A^ such that dip = (p. The polyhedron obtained is independent on the specific 
choice. 
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Figure 4. The dashed region represents the polyhedron V (Ji, 
when d{Ii,l2) = d > 0. In this case ^ "P (/i, I2). 
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2.3.3. Other topologies. We can generahze Theorem 12.141 to surfaces different from the 
two dimensional torus and to cellular decompositions different from the cubic one. More 
precisely consider an unoriented graph embedded on a compact surface without boundary. 
Vertices are associated to points of the surface, edges are associated to continuous self 
avoiding curves on the surface connecting vertices. Two different curves may intersect only 
on vertices. Cutting the surface along edges we obtain a finite number of two dimensional 
cells homeomorphic to two dimensional balls. Every edge belong to only two different two 
dimensional cells. From this we construct an oriented graph (V, E) whose vertices coincide 
with the vertices of the unoriented graph and whose oriented edges are obtained splitting 
any unoriented edge {x,y} into (x, y) and {y,x). Elementary cycles can be defined also 
in this case. They are naturally associated to one dimensional cells and oriented two 
dimensional cells and are defined like in the previous subsection. Also the boundary and 
co-boundary operators are defined in a similar way. We avoid formal definitions, see [3], 
[lOj and [14j for more details. We need to distinguish two cases, when the surface is 
orientable or not. If the surface is orientable and (j) E dA^ then ^ G A^ such that cf) = dtp 
is defined up to an additive constant. If the surface is non orientable and (/> G dA^ then 
ip & h? such that (f) = dijj Is uniquely determined. 

In the orientable case we can fix F'jq C F^ choosing an orientation for any two dimen- 
sional cell in such a way that any pair of adjacent elements of F'^ is oriented in agreement. 
We fix also E'^ C E^ choosing arbitrarily one among the two possible orientations for any 
edge. For any (x, y) G E'j^ there exists only one element /_|_ G F'j^ such that (x, y) G /+ 
and one element /_ G F'j^ such that (y, x) G /_. The corresponding interval I{x,y) is 
defined as in the case of the two dimensional torus. We then have the following Theorem. 

Theorem 2.15. Consider a weighted oriented graph constructed starting from a finite 
cellular subdivision of a compact orientable surface without boundary and let r £W such 
that (p^ = (p £ dh?' . Then we have r £ R^ if and only if 



Given Ji, . . . , J„ a collection of closed subsets of M we call V'{Ji, . . . , J„) C M" the 
closed unbounded convex polyhedron defined as follows. An element s = (si, . . . , Sn) G M" 
belongs to V{Ji, ..., Jn) if and only if the inequalities 

Si ^ d{0,Ji) , i = l,. . . ,n, 



16 D. GABRIELLI AND C. VALENTE 



are satisfied. Note that if Ji, . . . , J„ are closed intervals then the triangle inequality im- 
mediately implies V{Ji, ■ ■ ■ , Jn) ^ 'PiJi, ■ ■ ■ ■, Jn)- 

In the case of a non orientable surface we fix F'j^ C F/v choosing for any two dimensional 
cell arbitrarily one among the two possible orientations. Note that in this case is not 
possible to select the orientations in such a way that any pair of adjacent elements of F^ 
is oriented in agreement. Likewise we fix E'j^ C Ef^f choosing arbitrarily one among the 
two possible orientations for any 1 dimensional cell. 

For a non orientable surface, fixed {x,y) £ E'^, two possible situations are possible. 
The first case is when there exist one element /+ £ F'j^ such that (x, y) G /+ and one 
/_ G F'^ such that (y, x) G F'^. When this happens we define J{x,y) := I{x,y) as for 
orientable surfaces. The second case is when either do not exist elements of F'j^ to which 
(x,y) belongs or there are /i,/2 G F'j^j associated to adjacent two dimensional cells such 
that (x,y) £ fi, i = 1,2. We then define 

J{x,y) := \^-oo,mm {'ip{fi),^{f 2)} U max{'(/'(/i), V'(/2)} , +ooj , 

where depending on the cases either (y,x) G fi or (x,y) £ fi, i = 1,2. We have the 
following Theorem. 

Theorem 2.16. Consider a weighted oriented graph constructed starting from a finite 
cellular subdivision of a compact non- orientable surface without boundary and let r £ W 
such that (ff = 4> £ dh? . Then we have r £ R^ if and only if 

2.4. Applications. Finally we discuss some applications. In particular we concentrate on 
the elementary decomposition for the two dimensional torus. We start from a smooth con- 
tinuous vector field on T^ and consider its Hodge decomposition. The natural continuous 
counterpart of a discrete vector field on dh? is a vector field obtained as the orthogonal 
gradient of a smooth potential function tp. We introduce then a natural discretization 
procedure. The value of the discrete vector field on the edge (x, y) £ En is the value of 
the flux of the continuous vector field across the dual edge of (x, y) (see Remark [3] for 
the definition) with a normal vector oriented in agreement with (x,y). The result of the 
discretization procedure is an element of dA^ . Then we apply Theorem 12.141 to this very 
general example and obtain a condition for the validity of the elementary decomposition 
in terms of the variation of the potential function ip. 

Using the above framework we then construct a periodic random environment on Z . 
Theorem 12.141 gives a condition on the strength of the noise to be added in such a way 
that an elementary decomposition holds almost surely. Using the results in [5] we can 
then deduce a quenched Central Limit Theorem. 

All the applications are discussed in an informal way but the claims could be easily 
transformed into Theorems. 



3. Preliminary notions and results 

We start recalling some classic definitions of convex analysis. See for example [8j for 
more details. 
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Given w^, . . . ,w distinct elements of Z , with 2 ^ k ^ d + 1, we will say that they are 
in general position if the vectors 

/ r1 2 1 

a := w — w 

■■. (3.1) 

are linearly independent. By convention a single vector w"^ G IJ^ will be always considered 
to be in general position. It is easy to see that this definition does not depend on the 
specific order among the vectors. 

Given A C M*^ we denote by ^4'^ its relative interior part. This is defined as 

A° := {x E ^ : 3e > s.t. B,{x) n aff{A) C A] , 

where -Be(x) is the Euclidean ball of radius e centered at x. 

Given u; := {w^ , . . . ,w^) G {'Z'^)^ in general position then co{{w^, . . . ,w''}) is a (A; — 
1)— dimensional simplex and consequently for any y G co{{w^ , . . . ,w }) there exists a 
unique element fi := (^i, . . . , fik) of 

S'' := < ^= {fii,...,fik) ■■ l-H>0 , '^fii = l> , 



such that 



i=l 
k 



vO 



J2^^^w' = y. (3.2) 

i=l 

When y £ (^co({w^, . . . , w'^})) the coefficients /ij satisfy in addition the strict inequalities 
< /ij < 1. 

We recall that on W^, as shown in [16j . co{A) coincides with the set of elements x G M'^ 
that can be written as x = J^^ ydn{y), where ^ is any probability measure such that there 
exists a measurable set A' C A such that //(>!') = 1. In general this equivalence is false. 

Let us recall the following basic result of convex analysis (see for example 0) 

Theorem 3.1. [Caratheodory] Let A C M. . Then for any x G co{A) there exist 
x^, . . . ,x^ £ A in general position such that x G co{{x^, . . . , x^})^ . 

Remember that since the vectors are in general position then necessarily k ^ d+ 1. 
We state and prove the following characterization of the convex hull that we could not 
find in the literature. 

Lemma 3.2. Let x G M"^ and S C M'^. Then x G co{S) if and only if for any hyperplane 
H 3 X it holds 

H+nS^d) and H-nS^t/}. (3.3) 

Proof. Note that it holds {H^ Ci S)U {L[^ Ci S) = S. First we suppose that x G co{S) and 
show that (|3.3p holds. Assume by contradiction that for example there exists LI such that 
H^ n S = S and H^ n 5 = 0. This implies also that H~ n co{S) = 0. As a consequence 
we have 

X = X n co{S) CH-n co{S) = , 
a contradiction. 

Conversely we assume that ()3.3p holds for any hyperplane LL and show that this implies 
X G co{S). Assume by contradiction that x co{S). Then (see Theorem 4.4 in jBj) there 
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exists a separating hyperplane H among the two disjoint convex sets {j;} and co{S) for 
which it holds for example H^ 8 = 9 and x G H~^. If H is the hyperplane parallel to H 
and containing x then we have H~^ C H^ that implies H~^ n S = 0, a contradiction. D 

We recall the following basic result of convex analysis (see for example [8]) 

Theorem 3.3. [Helly] Consider {Ia}a^A ^ collection of compact convex subsets o/R"'. 
It holds Plae^ la 7^ 9 if end only if for any (ai, . . . , a^+i) G A'^^^ it holds ni=i ^Oi 7^ O- 

We defined the vector spaces A*^ , A^ , A^ associated to the discrete torus but it is clear 
that such vector spaces can be defined for any cellular complex. Correspondingly the 
following general result holds (see for example [3], [10] and [H 



Theorem 3.4. [Discrete Hodge Decomposition] For any finite cellular complex it 
holds the following orthogonal decomposition 

A^ = 5A° © dA^ © A\j , 

where A^ is called the subspace of harmonic one forms. 

We briefly discuss the ideas behind the proof of this Theorem considering the case of 
the 2 dimensional torus. The elements of 

(5A° := {</> G A^ : 3/ E A° s.t. = <5/} , 

are called potentials or gradient vector fields. The dimension of the vector space A*^ is 
|V/v|- It is easy to see that the kernel of the co-boundary operator 5 on A^ coincides with 
the constant functions, and in particular is a subspace of dimension 1. By the general 
identity 

dim (A°) = dim {Ker{S\^o)) + dim ((5A°) , (3.4) 

we deduce that dim{5h9) = \Vn\ — 1 = N'^ — 1, where dim{-) denotes the dimension. The 
orthogonal complement of ^A*^ in A^ is easily characterized. The subspace 5Js9 of A^ is 
spanned by the elements {51x}x<^v^ where Ix is the characteristic function of x G Vn- We 
deduce that an element (/> G A^ belongs to the orthogonal complement of 5h9 if and only 
if for any x G V/v it holds 

(0,5I,)i = -d</.(x) = O. 

This means that the orthogonal complement of 5Js9 is the set of divergence free discrete 
vector fields also called circulations. 

We have also that dim{h?') = \F'j^\ = N'^. It is easy to see that the kernel of the 
boundary operator d on A^ is the one dimensional subspace of the constant 2 forms 

Ker{d\^2) = {^ : V'(/) = c, V/ G F^ , c G M} . 

By the formula analogous to ()3.4p we deduce dim{d^) = N'^ — 1. To show a part of the 
orthogonal decomposition in Theorem 13.41 it is then enough to show that dip is divergence 
free for any t/; G A^. This is exactly the content of formula ()2.29p . Elements of dA^ are 
called 0— homologous circulations. The orthogonal complement (dA^) is characterized as 

7' 
N 



follows. Given g G Fj^ we define ipg G A^ as 



if / = 5 , 
^gif):={ -1 if/ = 5S 
otherwise . 
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Clearly dJS? is spanned by {dipg^gGF' • ^^ element (/> G A^ belongs to (dA^)-*- if and only 
if for any g S F'j^ it holds 

{cl,,d^g)i = 54>{g)=Q. (3.5) 

This means that the rotation of (f) around any g G F'j^ is zero. An element </> G A^ that 
satisfies condition ()3.5p for any g G F'j^ is called rotation free. Clearly any gradient vector 
field is rotation free. 
Finally if we define 

A]^ ■= (<5A°)^ n {dAY , (3.6) 

the orthogonal decomposition is proved. By (|3.6|) we have that the elements of A^ are 
rotation free circulations. Such a kind of discrete vector fields are called harmonic. Prom 
the dimensional counting we have dim (A^) = 2. 
Let us define (j)i £ A^, i = 1,2, as 

( +1 [iy = x + e^'^/N 
<l)i{x,y) := } -1 ify = x-e»/iV (3.7) 

[^ otherwise . 

As it is easy to check cpi and 02 are linearly independent rotation free circulations. Since 
dim (A^) = 2 we can identify 

A]j = {Ci0i + C2(j)2 , Q G M} . 

Consider a square matrix M whose rows and columns are labeled by a finite set V and 
call M{x, y) the element corresponding to row x and column y. The matrix M is called 
bi-stochastic if its elements are non negative and moreover 

l=Y,M{x,y')=J2M{x',y), yx,yeV. 

y'eV x'£V 

These conditions identify a convex compact subset called the Birkhoff polytope. To every 
element n G Sym,{V), the permutation group on V, we can associate a \V\ x \V\ matrix 
M,r called the permutation matrix. It is defined as 

It is clear that \V\ x \V\ matrices having positive elements and whose rows and columns 
are labeled by elements of V, are in bijection with weighted oriented graphs on V. The 
bijection identifies a matrix M and a graph {V,r) when M{x,y) = r{x,y) for any x,y £ 
V. For example the weighted graph corresponding to ()3.8p has weights r'^{x,y) := 1 if 
y = 7r(2;) and zero otherwise. The following result is classic (see for example [8]). 

Theorem 3.5. [Birkhoff- Von-Neumann] The set of\V\ x \V\ bi-stochastic matrices is 
convex and compact. Its extremal elements are the permutation matrices on V . 

4. Cyclic random walks on Z'^ 

Recall the definition of irreducible equivalence class of cycles stated just before Theorem 
12.41 The following Lemma identifies an important class of irreducible equivalence classes 
of cycles. 
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Lemma 4.1. Consider w = {w^, . . . ,w ) € (Z ) in general position and such that 
G co{{w^, . . . jw''})^. Then there exists an unique collection of strictly positive natu- 
ral numbers ni, . . . , rife such that 

[C] = {(u'^nl),...,('u;^nfc)} 

is an irreducible element ofC. We will call C* the set of irreducible elements of C obtained 
in this way. 

Proof. If y G co{{w'^ , . . . ,w^^)^ n 11^ then the element /i G S'^', uniquely determined by 
()3.2p has the coordinates jii that are rational numbers. Indeed if we call A the (d + 1) x /c 
matrix whose z-column is (lu^, . . . , w^, 1), then A has rank k. Every kx k sub-matrix 
containing the last row has determinant different from zero. This follows easily by the fact 
that the vectors are in general position. Moreover fj, is the unique solution of the linear 
system of d + 1 equations in the k variables //j 

k 

^ Aijfij = 6i^d+i + (1 - Si,d+i)yi > i = 1, . . . , d + 1 , (4.1) 

i=i 
where 5ij is the Kronecker delta. The unique solution to ()4.ip can be obtained by Cramer 
formula applied to k equations among which the last one. The matrices to be used have 
integer coefficients so that the solution is a vector of rational numbers. By construction 
they are also strictly positive. This means that if we fix y = in ()3.2p we obtain rational 
values for the //j. Let us write /ij = ^ where Oj and bi are natural numbers and ^ is an 
irreducible fraction for every i. Let b := lcm{bi, . . . , 6^} the least common multiple and 
define the natural numbers 



We have 



We show now that 



nj := 5/ii, i = 1,... ,k. 

k k 

^Uiw' = b^fiiw' = 0. (4.2) 

k 

Y^mw^^o, (4.3) 



j=l j=l 



when the strict inclusion 

{{w'^,mi),...,{w'',mk)} C {{w^,ni),...,{w'',nk)}, 
holds. Indeed, if ()4.3p is false then 

and since the vectors w^ are in general position we deduce 

^k = ^' = T- (^-^^ 



This implies 



lcm{bi, . . . , bk} ^ y^ mj < N^ Uj = b = lcm{bi, . . . , bk} , 



j=« j=« 
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and we get a contradiction. Equations ()4.2p and (|4.3p imply that if we define 

[C] :={(u;^nl),...,(u;^nfc)} , 

then [C] E C and moreover it is irreducible. 

It remains to show that the numbers Ui are uniquely characterized. Let us suppose 
that {{w^,mi), . . . , (w , m-fc)} is an irreducible element of C, then we want to show that 
necessarily rrii = rii. Clearly for such a collection of integer numbers rrii, equations ()4.4p 
and diS]) hold. We deduce Ej=i"ij = ^^> with / G N and / ^ 1. When I ^ 2 then 
{{w^,mi), . . . ,{w'^,'mk)} is not irreducible. When I = 1 then rrii = rii and we obtain 
[C]. D 

Not all irreducible equivalence classes of cycles belong to C*. Indeed consider 

«;i:=(_5,-5), w^:={l,l), w^ := (2,2) , 

vectors in 7? and the equivalence class {(w^, 1), (u;^, 1), (u;^, 2)|. It is easy to see that 
such equivalence class is irreducible nevertheless it does not belong to C* . 

Proof of Theorem ()2.3p . First we show that ii p £ B then it is cyclic. Since p £ B it holds 
()2.9p for any H. Equation ()2.9p can hold only if ()3.3p holds with S = S{p). By Lemma 
13.21 we deduce that G co{S{p)). We can then apply Caratheodory's Theorem 13.11 and 
deduce that there exists w = {w^, . . . ,w^) in general position with w^ G ^{p) ^iid such 
that G co{{w^ , . . . , w^})^ . Let [C] G C* be the corresponding irreducible equivalence 
class of cycles as constructed in Lemma 14.11 The corresponding purely cyclic measure is 

k 
i=l 

Here and hereafter we call /i(w) the unique element of S^ determined by (|3.2p with y = 
when the elements of w are in general position and moreover G co ({x^, . . . , x^}^ . 
Let 

mi:=mm\^^,i=l,...,k\>0. (4.7) 

We have that miqi :< p. This implies that pi := p — miqi belongs to TW ^ "^ and moreover, 
since qi is a mean zero probability measure, pi is balanced. Note also that by construction 
there exists x £ S{p) such that pi{x) = 0. Such an x is a vector w'^ that minimizes ()4.7p . 
We want now iterate this procedure. More precisely let pi be a balanced measure. 
Then G co{S{pi)). As before we can apply Caratheodory's Theorem identifying vectors 
w_ = {w^,...,w ) in general position, determine the corresponding element of C* and 
define like in (|4.6p the corresponding purely cyclic measure ^j+i. Then we define 

Pi{w^' 



rui+i := min <^ -—r- , j = l,...,k} > , (4.8) 

and finally call 

Pi+i := Pi - mi+iQi+i , (4.9) 

that is still a balanced element of 7W "^ ^. When |5(p)| < +oo, after a finite number / of 
iterations of the above procedure we obtain pi = 0. This follows directly by the fact that 
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l'5(Pi+i)l ^ I'5(Pj)I ~ 1- As a consequence we get 

I 
P = '^rniQi- (4.10) 

i=l 

Recalling that the probability measures qi are purely cyclic, equation ()4.10p is exactly a 
representation of p of the type 

p = pP= Y^ p([C7])p[^l, (4.11) 

[C]GC* 

where the measure p gives weight rrii to the unique element of C* associated to qi, i = 

When |5(p)| = +oo we need to implement the iterative procedure in a suitable way. 
Let us introduce the family of cubes 

A„ := < X G Z : max|xj| ^ n 

Given a balanced pi we define 

rii := inf {n e N : € co {S{pi) D A„)} . 

Clearly rij < +oo, since G co{S{pi)) and consequently there exists a finite number of 
elements in S{pi) whose convex envelope contains the origin. We can then define the 
measure g^j+i as in (j4.6p using vectors w. iii general position and belonging to A„^ n S{pi). 
Defining pj+i like in (j4.8p and ()4.9p we have 5(pj+i) C S{jpi) from which we deduce 
nj-|_i ^ n,. We obtain in this way an increasing family of cubes A„. C A„._|_^ and a 
decreasing family of measures Pi+i ^ Pi- We now show that necessarily 

hllli-J-+oo A„. = Z , (A'\')\ 

limj^+oo Pi = , 

where in the second limit it is enough to show just pointwise convergence. As a conse- 
quence, taking the limit j — t- +oo in 

i 
Pj=p-Y miQi , 

i=l 

we obtain 

+ 00 

P = ^rniqi. (4.13) 

i=l 

This identifies p with p^ like in (|4.1ip . where the probability measure p on C* gives weight 
rrii to the unique element of C* associated to qi. 

We need then to show ()4.12p . Both sequences are monotone and then the limits exist. 
Let us suppose by contradiction that limj_>.+oorai = n* < +oo. Note that in this case for 
any i we have S{qi) C A„* and moreover 

\S{pi+i)nAn*\ < \S{pi)nAn*\ -1. 

In particular after a finite number of iterations, say j we have that S{pj) D A„* = 0. This 
implies ?i-j+i > n* and we get a contradiction. 

Finally let us suppose by contradiction that liijii^^oo Pi = P* 7^ 0. Then, by mono- 
tone convergence Theorem, p* is balanced and consequently G co{S{p*)). Let us call 
w^, . . . ,w a finite collection of vectors in general position belonging to S{p*) and such that 



CYCLIC RANDOM WALKS 23 

S co{{w^, . . . , w'^})^. Consider j £ N such that {w^, . . . , w^} C An--i. By construction 
we have S{p*) C 5(pj) and moreover co (5(pj) n A„^_i). This is a contradiction. 

We now prove that if p is cychc then p £ B. By hypothesis (recaU ()2.8p we have 

p= hm y p([C])p[^], (4.14) 

[C]GC„ 

and for any fixed n the sum on the r.h.s. of (j4.14p satisfies (|2.9p for any H. We get 
immediately the vahdity of (|2.9p for p by applying the monotone convergence Theorem. D 

Proof of Theorem ()2.4p . The extremality of [C] E C* follows directly by the irreducibility. 
The validity of (|2.1ip has been implicitly proved during the proof of Theorem 12.31 D 

Proof of Theorem (j2.5p . To show the identification it is enough to show that TWq ^ 7W^ 
coincides with the set of cyclic measures. 

Let p be a cyclic probability measure, then clearly p £ A4^. Moreover it holds ()4.14p 

and this is the monotone limit assuring p £ A4q. 

Conversely take p G A^o H Ai^, we need to show that p is cyclic. By definition p G Ai^ 
and there exists a non decreasing sequence p" G A^o such that p = lim„^_|_oo p" . Since 
p" G TWo then also A" := p^ — p'^'^^ g A^Oj where we defined p^ := 0. This means that 
we can write A"' = ^^[(^ig^^* p"' {[C])p^'-^' for suitable measures /o". This follows by the fact 
that any element of A^o is balanced and then by Theorem 12.31 is cyclic. Since we have 
p = ^+~ A" we get dlSD with p{[C]) = ^+~ /o"([C]) and p on the l.h.s.. D 

5. Cyclic random walks on a finite set 

The following results are elementary but useful for the forthcoming results. 

Proof of Theorem. \ 2. 8\ Any purely cyclic graph is balanced. This implies that a necessary 
condition for the validity of (|2.14p is that r is balanced. Let us now show the other 
implication. Let 

m* := min r{x,y) > 

{x,y)eEir) 

and let (z^,z^) G E(r) such that r{z^,z^) = ?7i|. Due to the balancing condition 12.71 and 
the definition of tti^ there exists an edge {z^, z^) G E(r) with r{z^, z'^) ^ m^. By the same 
argument, if z^ 7^ z^, there exists {z'^,z^) G E{r) such that r(z^, z'^) ^ mj. We can iterate 
this procedure up to the first time we visit twice a vertex of V. Since V is finite this 
happens after at most \V\ iterations. We obtain in this way a sequence z^, z^, . . . , z"~^ of 
distinct elements of V such that (z*,z*^"^) G E{r) and moreover a z" such that z" = z^ 
for some ^ j ^ n — 1. We call Ci the cycle Ci := {z^,z^~^^, . . . ,z"). We also call 
nil := minj=j^. .. ^„_i r(z*, z*'^-'^) ^ m^ > 0. Clearly we have (V, m.iA'-"^') :< {V,r). As a 
consequence we have that r — miA'-^^' G W, it is still balanced and moreover it holds 
\E (r — mir'^^) | ^ E{r) — 1. This last inequality implies that after a finite number (at 
most l-E'(r)l) of iterations of the above procedure we obtain 

I 
r = y ruir^ *' , 

i=l 

that is the decomposition ()2.14p with the measure p that gives weight rrii to [Ci] £ C D 
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Remark 2. The validity of the balancing condition \2. 7\ implies that the uniform measure 
7r(x) := TyT satisfies the stationary condition 

7r(x) J2 r{x,y)= ^ 7r(y)r(y,x). 

{y:(x,y)eE{r)} {y:(y,x)eEir)} 

If the Markov chain is irreducible, this is the unique invariant measure. 

In the case of a continuous time Markov chain the set of balanced rates is convex but 
not compact. In the case of discrete time Markov chains, in addition to the balancing 
condition there are also the conditions ()2.12p . The constraints r ^ 0, ()2.12p and ()2.15p for 
every x, identify the Birkhoff polytope. The extremal elements of the Birkhoff polytope 
are then characterized by the Birkhoff- Von-Neumann Theorem 13.51 Theorem 13.51 together 
with the classical statement of Krein-Milmann Theorem [8j, ^^ implies that, given any 
bi-stochastic matrix M, we can decompose it like 

M = ^ m^M^ , 

TreSym{V) 

where m is a probability measure on SymiV). Written in terms of weights this equation 
becomes 



E 



ni-j^r'^ . (5.1) 

■jTGSym{V) 

Recall the classical result that every permutation can be decomposed into disjoint cycles. 
It is easy to see that in terms of weights this means that for every vr G Sym{V) we can 
write 

r-=J^r[^ri, (5.2) 

i 

where the C^ constitutes a family of disjoint cycles such that every element of V be- 
longs to one of them. The cycle containing the element x^ (^ V can be written as 
{x^ ,'k{x^),it'^{x^), . . . ,TT {x^),x^), where with vr*" we denote the composition of m— times 
the element vr € SymiV) and / is the minimal integer such that ^^^{x^) = x^ . Putting 
together (|5.ip and (|5.2p we obtain a special decomposition like ()2.14p . 

We finish the section discussing the case of an infinite graph. We simply reinterpret 
the results of subsection 12.11 in terms of an infinite weighted graphs with vertices T/ = Z 
and edges i? = Z x Z . The weighted graph corresponding to the translation invariant 
Markov chain on V = U^ defined by ()2.3p . gives weight r{x,y) = p{y — x) to the edge 
{x,y). Let us suppose that p is a cyclic measure. Consider Cj a cycle representant of the 
equivalence class in C* corresponding to the cyclic measure qi in ()4.13p . On Z there is 
defined a shift operator Tx that acts naturally on cycles by 

I ^\tiy % Jb ^•••^>Zj ^ Jb 1. I Jb I Jb t^ Jb \ Jb ^ , , , ^ Jb \ Jb ^ Jb \ Jb J • 

Define the family of cycles {TxCiY^^-^d in U^. Let p be the positive measure on C that 
gives weight rrii to [r^^Cj] for any x. The results in subsection 12.11 implv the validity of the 
decomposition ()2.14p with this specific measure p. In this case ()2.14p becomes 

r(x, y) = p{y - x) = ^^ niiA'^''^'^ {x, y) , V (x, y) . 
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6. Cyclic random walks on a finite graph with topology 

We discuss here the general case of a d ^ 2 dimensional torus. More detailed results for 
the specific cases d = 1,2 will be discussed separately in some subsections. 
We start observing that formula (|2.25p can be written as 

A subset A C VF is monotone non decreasing if r £ A and r < r' implies r' £ A. The 
subsets R^ and R* are in general not monotone subsets of W . Nevertheless R^r]R{(j)) and 
R* n i?(</>) are non decreasing for any fixed (j). This is the content of the next lemma. 

Lemma 6.1. Consider r, and r' belonging to W and such that r ^ r' and (jf = (jf . We 
have that if r £ R* then also r' £ R* . The same happens for R^. 

Proof. We prove the statement for R* . The proof for R^ is the same. By assumption there 
exists a decomposition like (I2.21|) for r. We then have 

r' = Y, p([C])r[^]+(r'-r). (6.1) 

[C]ec- 

By the hypotheses of the lemma h := r' — r belongs to W and moreover (p = 0. Clearly 
R{0) ^ R* and moreover for any h £ R{0) we have 

h= Y^ h{x,y)A^i-'y}^ . (6.2) 

Putting together ([61]) and ([O]) we get r' £ R*. D 



Proof of Lemma \2.y\ Since we discussed the discrete Hodge decomposition in the two 
dimensional case we will prove this Lemma also in the two dimensional case. The proof 
in the general case is analogous. Consider C = {z^, z^, . . . , z") with z"^ = z^ a, cycle such 
that [C] £ C* . Since C is homotopically trivial it holds (j2.20p . Componentwise ()2.20p is 
written as 

|{(x,a; + eW/iV) G C7}| - |{(x + e«/iV,2;) gC7}| =0, 
|{(x,x + e(2)/Ar) £ C}\ - \{{x + e^^^ / N , x) £ C}\ = , 



that implies 






where the vector fields (pi are defined in (j3.7p . This means that (j)^ £ (A^) . More- 
over we know that A'-"' satisfies the balancing condition 12.71 This implies that 4>^ is a 
divergence free discrete vector field i.e. it is a circulation. This is equivalent to say that it 
belongs to ((5A°)^. By Theorem lOl we obtain </>'''''' £ dK^ for any [C] £ C* . Given r £ R* 
having a decomposition like ()2.2ip . by linearity of the projection of weights onto vector 
fields we have 

Since dA^ is a vector subspace, it is closed under linear combinations and consequently 
the right hand side of ()6.3p belongs to dA^. Since C^ C C* the condition (j)^ £ dh?' is also 
necessary for the validity of (|2.22p . D 
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Figure 5. An r £ W for the two dimensional torus such that 0'" G dA^ 
nevertheless r ^ R*. Positive unitary weights are represented by boldfaced 
arrows. Opposite sides of the square are identified. 



Proof of Lemma \2.1U\ We prove the statement for R* . The proof for i?^ is the same. 
We need to prove that for any sequence r„ E R* converging to some element r £ W we 
have necessarily r £ R*. Observe that for any {x,y) G Ejy we have rn{x,y) ^ and 
consequently linin^^oo rn{x , y) = r{x,y) ^ 0. Since for every n it holds r„ G R* , we can 
write 

for some positive measures pn- Let M := maxj-^. y)^^;^^ r{x,y). Consider any [C] G C* and 
take (x, y) G C. For n large enough we deduce 

O^PniiC])^ Y, Pn([C])r[^](x,y) = r„(x,2/)^2M. 
ic]ec'- 

The last inequality follows by the fact that r„ converges to r and r{x,y) ^ M. Since 
Pn{[C]) takes values on a compact set we can then extract a converging subsequence 
PrijiiC]). By a finite Cantor diagonalizing argument, there exists a subsequence (that we 
still call rij) such that Pnj{[C]) is converging for any [C] G C*. Let us call /o([C]) ^ the 
corresponding limits. Since we have a finite sum we get 

r = lim Vn, = lim ^^ «-- r[r71V^ I 



lim 



pn.my 



[C]GC* 



E 

[C]ec' 



p([C])r 



[C] 



D 
" G dA^ 
but the 



This is equivalent to say r £ R*. 

The bi-dimensional example illustrated in figure [5] shows that the condition ^ 
is not sufficient neither for ()2.2ip nor for ()2.22p . Indeed in this case (p"^ £ dh? 
oriented weighted graph (Vat,?') contains no homotopically trivial cycles. Consequently 
decompositions of the type ()2.2ip and (|2.22p are not possible. 

Recalling Lemma F6.1l and the fact that r'*' is the minimal element of R{4')i it is natural to 
study the following problem. Given (p £ dl\^ a discrete vector field, under which conditions 
r't> £ R*, or r"^ G R''? When r<^ G R* then by Lemma O we deduce that R{(p) C R* . 
When r^ £ R^ then by Lemma l6. II we deduce that R{4') C. R^ <^ R*. On the other side 
if r"^ R* then by Lemma 12.101 all the elements of VF in a neighborhood of r"^ will also 
not belong to R*. Note that inside such a neighborhood there will be also rates such that 
r{x,y) > for any {x,y) £ E^. The same happens for R'^. 
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Proof of Corollary \2.11\ The proof follows directly from the above argument and the fact 
that there exists a (^ G dA^ such that r'*' ^ R* . In dimension two this is exactly the example 
in Figure O Similar examples can be clearly constructed in any dimension d ^ 2. D 

We will give a complete characterization of vector fields such that r*^ G R^ in dimension 
one and two. A characterization of vector fields such that r*^ G R* seems to be an 
interesting combinatorial problem. 

6.1. One dimensional torus. The balancing condition (j2.15p at site x G Vn can be 
written as 

r{x, X + 1/A^) — r[x + l/N, x) = r{x — 1/N, x) — r{x, x — l/N) , 

and corresponds to require that the discrete vector field (p"^ has divergence zero. Note that 
in one dimension a discrete vector field is divergence free if and only if it is constant. 

Proof of Theorem \2.12\ We want to characterize the Markov model having rates r admitting 
a decomposition of the type 

r = E P(iC])r^''^ = Y. P([C{x,x+im])^"'^^'^+^^"*' • (6-4) 

If ()6.4p holds then clearly we have that the associated discrete vector field satisfies 

(/>''(X,X + 1/N) = p([C|^.,a.+l/Af}]) - P{[C{x,x+l/N}]) =0. 

Conversely given a weight r such that cj)'^ = then we can set in ()6.4p p{[C^x^x+i/N}]) = 
r(x,x + 1/A^) = r{x + 1/N,x). Finally the reversibility follows by the detailed balance 
condition 

—r{x,x + l/N) = —r{x + l/N,x), Vx £ Vn , 

that is equivalent to (j)^' = 0. D 

Proof of Theorem \2.1S\ We want to characterize now the Markov models having a decom- 
position of the type 

r= E ^([C])^'^^= E /^([q.,.+i/iv}])r[^^--+^/->l+p([C+])r[^+l+p([C_])r[^-l. (6.5) 
IC]£C xGVm 

The discrete one dimensional torus is a finite graph and consequently this problem has 
been solved in Theorem 12.81 We discuss it again in terms of the discrete vector field <j)'^ . 
The fact that in one dimension a discrete vector field has zero divergence if and only if it 
is constant has already been stressed. If (|6.5p holds then clearly we have 

0^(x,x+i/iv) = p([q,,,_,i/^}])-p([q,,,+i/^}])+p([c+])-p([c_]) = p([c+])-p([c_]), 

that is a constant discrete vector field. Conversely consider some rates such that for any 
X G V/v it holds (f)^{x,x + 1/A^) = c with c a constant real number. Then from ()6.5p 
recalling the decomposition ()2.25p we obtain 

r(x, X + 1/N) = p{[C{,^x+i/N}]) + p{[C+]) = [c]+ + s{x, x + 1/N) , 

r(x + 1/N, x) = p{[C{,^x+i/N}]) + pi[C-]) = [-c]+ + s{x + 1/N, x) , ^^'^^ 

where we recall that s is the symmetric part of r. Clearly (|6.6p implies 

p{[C{x,x+i/N}\) = s{x, X + 1/N) - a , (6.7) 

where a is a real parameter. Putting ()6.7p in (|6.6p we obtain also 

p{[C+\) = [c]+ + a, , 
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Since the left hand sides of ()6.7p and ()6.8p are non-negative we obtain the constraint 
o G [0, 77i]. Recahing the definition of the symmetric part s we obtain ()2.30p . D 

Note that the decomposition ()6.5p is unique only when tti = 0. 

6.2. Two dimensional torus. Let p a positive measure on C^. The general decomposi- 
tion (El2\\ reads 

r = E P(iC{.,y}])r^''^^'^'^ + E PiiCf])r^'''^ ■ (6.9) 

To the measure p in ()6.9p we can also associate an element tp^ S A^ defined as 

rif) ■■= p{[Cf]) - p{[Cf.]) feFM. (6.10) 

Lemma 6.2. If r is like in (j6.9p then we have (j/' = dtp^ . 

Proof. Given {x,y) G Ejsi it belongs to only two elements of -F/v- Moreover one is clockwise 
oriented and the other one is anticlockwise oriented. Let us call /+ the element of F'^ 
such that (x,y) £ /+ and /_ the element of F'^ such that {x,y) G fl. Then we have 

f{x, y) = (p([C/j) - p([C/c ])) - [p{[Cf_]) - p{[Cr_])) = nu) - rif-) . 

Recalling ()2.27p and ()6.10p we obtain the statement of the lemma. D 

The following Remarks will not be used during our proofs but will be useful in the 
following. 

Remark 3. (Duality) We recall a well known duality relationship. To the 2 dimensional 
discrete torus (YnjEn) we associate a dual discrete torus (VntE^) defined as follows. 
The vertices of the dual graph are the elements ofT^ having coordinates x-\-2j^{e^^' +e'^^) 
with X G V/v- Note that every element of Vjv is the center of a cell of the original cellular 
decomposition ofT"^. The set of oriented edges E^ is constituted by the pairs {v,w) such 
that v,w £ Vj\f and moreover d{v,w) = 1/N. It is possible to define a duality map D. 
This map is defined both on E^ with image on E^ and on E^ with image in Ej^j. This 
map is infective and satisfies the involution property D^ = —I. We use the same symbol 
D, both when it acts on Ej^ or in E^, since it can be defined easily in the same way. Any 
element {x,y), both of E^ and of E^ can be naturally represented by an arrow exiting 
from X and entering in y. The element D{x,y) is defined as the unique element of either 
En or En, whose representing arrow is obtained rotating counterclockwise of ^ the arrow 
representing (x, y) around its middle point. The dual map can be naturally extended to act 
on discrete vector fields. More precisely given for example cj) £ A^ we define D<p G A^ as 

{D(j)){w,z) := (l){x,y), 

where {x,y) is the unique element of E^ such that D{x,y) = {w,z). Clearly we called 
A^ the set of discrete vector fields on the dual torus. According to these definition the 
notion of circulation and rotation free are dual to each other. More precisely given a 
discrete vector field 4> whose discrete divergence is zero in a vertex, then D(j) will have zero 
circulation around the corresponding dual face. Conversely given a discrete vector field (p 
that satisfies the condition of zero rotation around a face then Dcf) will satisfy the condition 
of zero discrete divergence in the corresponding dual vertex. 
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Remark 4. Given tp ^ A'^ it is simple to compute (j) = dip by 

(t>{x,y)=i;{U)-ij{f.). (6.11) 

In terms of the dual graph it can be interpreted in the following way. Consider ip the 
element of Aq that associates to any vertex of Vn the value ip{f) where f G Fj^ is the 
oriented dual face of the fixed vertex. Then (|6.1ip can be written as 

D(/) = 6'tp. (6.12) 

In extended form (|6.12p is 

{x, X + e(i)/7V) = dipix + e(i)/(2iV) - e(2)/(2iV), x + e(i)/(2Ar) + e(2)/(27V)) 



p{x, X + e(2)/Ar) = -d4,{x - e^^^/{2N) + e^^^/{2N),x + e^^^/{2N) + e(2)/(2iV)) . 

(6.13) 
Conversely given cj) G dA^, the determination of a ip such that (p = dip requires a non local 
computation. More precisely using (|6.12|) we determine 5ip G Ai and then we can compute 

n 

iP{f)-iP{f') = Y,5ij{x\x'+^), (6.14) 



4 = 1 



where (x*,x*+ ) G E^ and x is dual to f and x" is dual to f. 

Proof of Theorem \2.14\ Given cj) G dA^ there exists a ip G A"^ such that (p = dip. Since 
the kernel of the boundary operator d coincides with the constant 2 forms, the elements 
ip' G A^ such that dip' = cp are exactly of the type ip + c where c is an arbitrary constant. 
Let r £ W having a decomposition like (j6.9p and such that 0'" = (p. By Lemma (|6.2p 
we have ipP = ip -\- c for some constant c. This means 

p{[Cf]) - p{[Cf^]) = iPif) + c , / G F^ . (6.15) 

Recalling ()2.25p we have that 

s:=r-r'l' (6.16) 

is the symmetric part of r and belongs to i?(0). We want to characterize which are the 
elements s G i?(0) that can be obtained in ()6.16p when r G i?*^ n R{(p). 
The more general positive solution to (|6.15p is 

P{[Cf]) = [{ip + c)(/)]+ + m{f) , feF^, (6.17) 

where m{f) are arbitrary non-negative numbers such that m{f) = m{f^). In (j6.17p ip + c 
is the element of A^ obtained by the sum of iJj and the constant 2-chain c G A^ defined by 
c(/) = c for any / G Fj^. In particular (j6.17p means 

PiiCf]) = mf) + c]+ + m(/) , p{[Cf.]) = i-iPif) - c]+ + m(/) , / G F^ . 

Note that we can write 

Y^ m(/)r[^/] = J2 \ Yl "^(-^^ I ^'^*"'"*' • (6-18) 

Putting (I6T7)1 in (^ using ([638]) we get 

« = E /''([qx,,}])r-[^{-^>] + E [(V' + c)(/)]+^t'^^l - r'f' , (6.19) 
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where 

p'([q.,y}]:=p([q.,y}])+ E "^(/)- 

{f€FN ■■ ix,y)€f} 

Since the positive numbers m{f) and p{[Cx,y]) are arbitrary, for any fixed c E M (|6.19p 
says that s ^ s{c) where 

s(c):= J^[(V + c)(/)]+r[^/]-r^ 

Recall that given (x,y) G £'^ there exist only two elements of Fjy that contain {x,y). 
One is clockwise oriented and the other one is anticlockwise oriented. As before we call /+ 
the element of F^ such that {x,y) G /+ and /_ the element of F^ such that (x,y) G fi. 
If c G M is fixed, the condition s ^ s(c) is equivalent to require for any (x, y) G E'j^ (the 
condition for (y, x) will automatically be satisfied) 

s{x, y) ^ [iPiU) + c]+ + [-V(/-) - c]+ - [<t){x, y)]+ . (6.20) 

Since (p = d^ we have also 

0(x,y)=V(/+)-V'(/-). (6.21) 

It holds the following identity 

where a and 6 are real numbers, sg(-) is the sign function that associate to any real number 
its sign and the first identity is the definition of the function Z. Putting (|6.2ip in (|6.20p 
we get for any (x, y) G E'j^ 

s{x,y)^Z{i;{U) + c,i^if.) + c). 

It is simple to check that 

Z{i;{U) + c, V(/-) + c) = d(0, c + I{x, y)) , 

where d is the Euclidean distance on the real line and I{x, y) is the closed interval of the 
real line [min{V'(/+), V'(/-)} ,max{^(/+), V'(/-)}] := [ii{x,y),i2{x,y)]. 
Let us call 

S{c) := {s G R{0) : s ^ s{c)} . 
We showed that 

S{c) = {se R{0) : s{x, y) ^ d{0, c + I{x, y)) , {x, y) G E'^] . 

Since c is arbitrary we get that when r G M*^ n i?(</>) the corresponding symmetric part s 
in ()6.16p satisfies 

s G UegR5(c). (6.22) 

We now write this condition in a simpler form. Indeed we will now show that s G -R(O) 
belongs to UcgR'S'(c) if and only if 

^(x,y)eE'^ [-i'2{x, y) - s{x, y), -h{x, y) + s{x, y)] / . (6.23) 

Indeed it holds 

s(x,y) ^ d(0,c + /(x,2/)) ^^ c G [-i2(x,y) - s(x,y), -ii(x,y) + s(x,y)] . 

(6.24) 
If s G Uc5'(c) then, using (|6.24p . there exists a c* G M such that 

c* G [-i2{x, y) - s(x, y), -ii{x, y) + s{x, y)] , V(x, y) £ E'^ , 
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and ()6.23p holds. Conversely assume that ()6.23p holds and call c* an element of the non 
empty intersection. Then still using ()6.24p we get s G S{c*). 

Since closed intervals are compact convex sets of M^, we can write condition (|6.23p in a 
simpler form using the Helly's Theorem 13.31 We apply it in the special case of dimension 
one. We obtain that condition ()6.23p holds if and only if for any pair {x,y) and {x',y') of 
elements of E'j^ it holds 

[-i2{x,y) - s{x,y),-ii{x,y) + s{x,y)] D [-i2{x',y') - s{x' ,y'),-ii{x',y') + s{x',y')] / 0. 

This is equivalent to require 

d{l{x,y),I{x',y')) ^s{x,y) + s{x',y'). (6.25) 

Since ()6.25p has to be satisfied for any pair of elements of E'j^ we deduce that condition 
(|6.22p is equivalent to 

D 

Let us briefly discuss how the above results can be used on a problem concerning an 
infinite graph. Consider the infinite graph having vertices Z^ and edges coinciding with 
the pairs of vertices {x,y) such that d{x,y) = 1. Consider also a periodic weight r, i.e. a 
weight such that there exist A^i and N2 integer numbers such that 

r{x, y) = r{x + (iVi, iVa), y + {Ni,N2)) , V(x, y) . 

Then the problem of finding a periodic cyclic decomposition of r, i.e. a decomposition such 
that /o([C]) = p([T(Afi,Ar2)C'])) is strictly related to the problem of finding a decomposition 
like ()2.2ip or ()2.22p for a finite torus with Ni x N2 vertices. Indeed in the case of elementary 
decompositions ()2.22p the two problems are equivalent. Apart from an irrelevant dilation, 
an elementary cycle on the discrete torus can be naturally interpreted as an elementary 
cycle on the infinite graph. It is easy to see that if ()2.22p holds on the discrete torus, then 

E E ^([^])^'"''^ (6.26) 

is a periodic cyclic decomposition on the infinite graph with elementary cycles. In (|6.26p 
C^ are the equivalence classes of elementary cycles on the discrete torus. Conversely every 
elementary cycle on the infinite graph can be projected onto an elementary cycle of the 
torus. By periodicity of the weight r you get a decomposition like (j2.22p . In the case of 
decompositions ()2.2ip the relation among the two problems is more subtle and we will not 
discuss here. 

6.3. Other topologies. The proofs of Theorems 12.151 and 12.161 are very similar to the 
one of l2.141 We give only a sketch. 

Proof of Theorem \2.15\ The proof of this Theorem follows closely the lines of reasoning of 
the proof of Theorem 12.141 In particular the two basic ingredients are the following. The 
first one is that every edge {x,y) E E'at belongs to only two elements of F/y. The second 
one is that we can choose orientations in F'^ in such a way that every pair of elements on 
-F^ is oriented in agreement. As a consequence given (j) S dl^ the two chain ip such that 
(j) = dip is defined only up an additive constant. D 

Proof of Theorem \2.16\ The proof of this Theorem uses the same type of arguments of 
Theorem 12.141 The difference is that given (p G dA^ then the two chain ip such that 
(p = dip is uniquely determined. This is due to the fact that the surface is not orientable. 
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For any (x, y) £ E'j^ such that there exist f+, f- G Fj^ such that (x, y) G /+ and (y, x) G /_ 
the we get the constraint s{x,y) ^ d{0,I{x,y)). For an {x,y) £ E'^ for which for example 
there exist /i, /2 G i^^ such that (x, y) £ fi, i = 1, 2, instead of (|6.20p and (|6.2ip we get 

six, y) ^ [^(/i)]+ + [V(/2)]+ - [^{h) + V'(/2)]+ . 
Since for any pair of real numbers a, b it holds 

H+ + [6]+ - [a + 6]+ = I --il«l' 1^1} if ^g(«) ^ «g(J) ' 
^ ^^ ^ ^^ ^ ^^ [ if sg(a) = sg(o) , 

we deduce in this case s{x,y) ^ d{0, J{x,y)). The Theorem follows. D 

7. Applications 

We start this section discussing a very general example. First we show how to discretize 
in a natural way a smooth divergence free vector field by a divergence free discrete vector 
field. Then we apply the results of Theorem 12.141 In dimension d = 2, the continuous 
version of the Hodge decomposition is the following. Let u £ T"^ and F = (Fi{u),F2{u)) be 
a smooth vector field on T^. The Hodge decomposition says that there exist two smooth 
functions / and Tp such that 

F = V/ + V^V + ai(l,0) + 02(0,1). (7.1) 

In the above formula V is the gradient, V^ip{u) := {tlju2{u), —ipui{u)) is the orthogonal 
gradient and aj = Jj.2 Fi{v)dv. The definition of the orthogonal gradient can be seen as 
a continuum version of formula (|6.13p in Remark [H Decomposition (|7.ip says that the 
continuous version of the discrete vector fields on dA^ are the vector fields on T^ of the 
form 

F={iPu„-4^u,), (7.2) 

for a suitable smooth function ip. Clearly a vector field of the type ()7.2p has zero divergence 
since 

V • V^ip = ipu2Ui - ipuiU2 = . (7.3) 

To any smooth vector field of the type ()7.2p we can associate a discrete vector field (pN £ 
dA^ on the discrete torus defined as 

(/)N{x,x + e^ '/N) := ^ ^U2 {xi + —,y j dy 

= V'(xi + ^,X2 + ^)-^(xi + ^,X2-^), (7.4) 



and 



0Af(x,x + e(2V^) := - / ^ i^u^iy,X2 + —]dy 

"2N ^ 



X\- 



= V'(^l-^,X2 + ^)-^(xi + ^,X2 + ^). (7.5) 

In ()7.4p and ()7.5p we fixed the values of the discrete vector field coinciding with the fluxes 
of the continuous vector fleld V V respectively across the segments 



c 



[x + e^^^/{2N) - e(2)/(2iV)) + (1 - c) (x + e^^y{2N) + e(2)/(2iV)) , c E [0, 1] , 
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and 

c (x - e^^'^/{2N) + e(2V(2A^)) + (1 - c) (x + e^^'^/{2N) + e(2)/(2iV)) , c e [0, 1] , 

with associated normal vectors respectively e^^' and e^"^'. According to these definitions 
it is easy to find a two chain if^^ G A^ such that (f)^ = dipN. Given / G F^ of the form 

/ = {x,x + e^'^'^ /N,x + e^^^ /N + e^'^'^ /N,x + e^'^'^ /N,x) we have 

The condition of zero divergence for (j)^ can be also directly derived by (|7.3|) . 
We us call 

M := max \ip{u) — il^{v)\ . 

Fix a weight r := r'^'^ + s G R{4>n) with s{x,y) ^ — for any {x,y) G En. Since 
d{I{x,y),I{x' ,y')) ^ M for any pair of elements of E'j^ then 

s{x, y) + s{x', y') ^ M ^ d{I{x, y), I{x',y')) , 

so that r e R^. 

As an application of this example we discuss a direct consequence for a random walk 
in a random environment. Consider the infinite graph with vertices I? and oriented edges 
coinciding with pairs of nearest neighbor vertices. Let ^ be a smooth periodic function on 
M^ with integer periods A'^i and N2, i.e. such that ipiu) = ip{u+ {Ni, N2)) for any it G M^. 
Let U = {Ui,U2) be a random variable uniform on [0, A^i] x [0, A'^2] ^ 1^^- We construct 
the random discrete vector field (j){uj) defined by 



(x,x + e(^);w) := Tu(uj) 



(x,x + e^'^';uj 



ru{u) 



V'l ^1 + 2'^^"^ 2) ~^ (^^"^ 2'^^~ 2 

W [Xi- -,X2 + -] - 1p [Xi + -,X2 + 



2 2J ^ \ 2 2, 

Clearly 4>{u:) G dh?' for any cj and moreover the law of (f) is invariant under translations. 
Let also {U{x,y}\i \^^i be a collection of i.i.d. random variables taking values on an 

interval [a, 6], a, 6 > and independent from U . Here E'^ are oriented edges of a A^i x A'^2 
discrete torus of mesh size A^ = 1. We construct the random weights 

r(x,y;w) := Z{x,uj)-^ (r'^('^)(x,y) + C/|,,,|) , (7.6) 

where Z[x,uj) is the normalization constant that guarantees the validity of ()2.12p . Note 
that C^{x,y} = U^x+Ni,y+N2} so that r{oj) is a periodic enironment. Let 

M := max \ij){u) — 4'{v)\ ■ 

u,v<E[0,Ni]x[0,N2\ 

If a ^ — then the hypotheses in [5] are satisfied and the discrete time random walk 
in random environment having transition probabilities determined by (j7.6p satisfies a 
quenched Central Limit Theorem since it admits a periodic cyclic decomposition into 
elementary cycles. Since the random environment is periodic the result obtained is rather 
simple but the argument can be generalized to the non periodic case [7|. 

We end the paper with some comments. When r G W^ is given, to determine if r G i?^ 
we need to compute the associated discrete vector field (j) and then a two chain ip such that 
(j) = dtp. Once ip is computed we determine the corresponding polyhedron V and finally 
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we can check if the symmetric part s of r belongs or not to V. As discussed in Remark U] 
the determination of ip requires a non local computation. It is then useful to have some 
simple sufficient conditions. For any pair {x,y) and {x',y') of elements of E'^ we have 

d{I{x, y),I{x',y')) ^ max \i;{f) - V'(/')| • (7.7) 

If we are able to estimate from above the right hand side of (|7.7p by a constant M and 
s{x,y) ^ M/2 for any {x,y) £ E'j^ then r G R'^. Consider ip the element of Aq associated 
to ip as in Remark H] and use (|6.14p . Consider an unoriented graph having as vertices Va^ 
and as unoriented edges {x,y} with {x,y) G E'jyj. To every unoriented edge we associate 
the weight w{x,y) := \(jf{D{x,y))\. The weight of a path is the sum of the weights of 
its edges. These weights introduce a metric structure on the graph. The distance among 
two vertices is the minimal weight among all paths connecting the two vertices. An upper 
bound to the right hand side of (|7.7p is then the diameter of this unoriented graph. It can 
be estimated for example as the weight of any spanning tree. 
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